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We consider an arbitrary continuous variable three-party Gaussian quantum state which is used 
to perform quantum teleportation of a pure Gaussian state between two of the parties (Alice and 
Bob). In turn, the third party (Charlie) can condition the process by means of local operations and 
classical communication. We find the best measurement that Charlie can implement on his own 
mode that preserves the Gaussian character of the three-mode state and optimizes the teleportation 
fidelity between Alice and Bob. 
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I. INTRODUCTION 



Quantum information processing with continuous variables (CV) provides an interesting alternative to the tra- 
ditional qubit-based approach. CV seem to be particularly suitable for quantum communications, as for example 
quantum teleportation pj and quantum key distribution 0. Multipartite CV entangled states for quantum commu- 
nication networks are rather easy to produce. In particular, tripartite entangled Gaussian states realized either using 
squeezers and beam splitters interlinked bilinear interactions Q| or through radiation pressure HUES 13 have been 
considered for the realization of a quantum teleportation network and for telecloning. Tripartite CV entanglement 
between optical modes has been generated using squeezing and beam splitters and experimentally characterized in [j| , 
and it has been recently exploited for the realization of quantum secret sharing in and for quantum telecloning 
in • Here we consider a generic CV tripartite Gaussian state which is employed for the specific task of teleporting 
a pure Gaussian state between two of the three parties (Alice and Bob) . We determine the best way the third party 
(Charlie) can cooperate to improve this teleportation task. To be more specific, we find the optimal Gaussian mea- 
surement at Charlie's site which maximizes the teleportation fidelity. This is different from the optimization over all 
possible local Gaussian operation of CV teleportation as considered in ^1 an d a l so from the problem of entanglement 
distillation |l3j . where one always starts from bipartite entangled states and tries to increase their entanglement. In 
Sec. II we present the scenario and describe the teleportation protocols in the case when it is assisted or is not assisted 
by measurements at Charlie's site. In Sec. Ill we discuss the case when Charlie performs a dichotomic measurement 
with a Gaussian and a non-Gaussian outcome, while in Sec. IV we consider the case of a local Gaussian measurement 
performed at Charlie's site. In Sec. V various applications of the theorems derived in Sec. Ill and IV are discussed 
in detail, while Sec. VI is for concluding remarks. 



II. ASSISTED AND NON-ASSISTED TELEPORTATION PROTOCOLS 

The scheme we are going to study is described in Fig. 1: Alice, Bob and Charlie each possess a continuous variable 
mode, characterized by an annihilation operator a, b and c respectively, and share a quantum channel given by an 
arbitrary three-mode Gaussian state p, characterized by a displacement dgK 6 and a correlation matrix (CM) 

/ A F E\ 

V = \F T B D , (1) 
\E T D T C ] 

where the blocks A,B, ...,F are 2x2 real matrices. Alice has to tclcport to Bob an unknown pure Gaussian state 
Pi n with CM Vi n and amplitude p. The most straightforward strategy is to ignore Charlie (non-assisted protocol) 
and use the reduced bipartite Gaussian state p tr = Tr c (p) to implement a standard continuous variable teleportation 
protocol 0. In such a case, Alice mixes her part of the reduced state p tr with the input state pi n through a 
balanced beam-splitter and makes a homodyne detection of the output modes, i.e., she measures the quadratures 
X- = 2~ 1 / 2 {X a — Xin) and P+ = 2~ 1 / 2 (P a + Pi n ). After the measurement, Alice classically communicates the result 
7 = — X- + iP + to Bob, who performs a conditional displacement Dbil') = exp(7'6^ — j'*b) on his own mode 6, 
where 7' = 7 + S has the double effect to compensate the displacement due to Alice's measurement (by 7) and the 
displacement of the reduced state (by S which is connected with d, see ) ■ This means that Bob can always implement 
(through a suitable displacement 7') a displacement-independent teleportation protocol, whose fidelity only depends 
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on the CMs of the reduced state and the input state. When Charlie is traced out, p tr has the CM 

and the teleportation fidelity is given by F tr = (detr^)" 1 / 2 H2, where 

T* r = 2V m + RAR + B-RF- F T R (3) 

and R — diag (1,-1). 

An alternative strategy for Alice and Bob is to ask for the help of Charlie (assisted protocol), who can perform a 
suitable measurement on his own mode c and classically communicate the result to Bob (see Fig.0. In this modified 
protocol, Bob performs his displacement only after receiving the information about the measurement outcomes from 
both Alice and Charlie. For every Charlie's outcome n (with probability P n ), Bob can choose a conditional displacement 
7^ aiming at optimizing the conditional fidelity F^ n > and therefore the effective fidelity F = ^2 n P n F^ nS> of the protocol 
[l4j . In particular, if the bipartite reduced state conditioned to the outcome n, p^ n \ is a Gaussian state, then Bob's 
displacement is given by j' n = 7 + S n where S n exactly cancels the displacement of p^ , and therefore the conditional 
fidelity F^ n > depends, as before, only on the CMs. In the following we consider two general kinds of measurement 
at Charlie's site: a local dichotomic measurement, with a Gaussian outcome and a non-Gaussian one, and a local 
Gaussian measurement, defined as a local measurement preserving the Gaussian character of the shared state for 
every outcome. Our aim is to compare the assisted fidelity F and the non-assisted fidelity F tr for both kinds of 
measurement. We anticipate that for the dichotomic measurement one does not have an improvement (F < F tr ), but 
the results achieved for the conditional fidelities F^ are interesting and they can be directly extended to the case of 
the Gaussian measurement where one can optimize F and then surely state that F > F tr . 




FIG. 1: An arbitrary 3-mode Gaussian state is shared by Alice, Bob and Charlie. Alice is supplied with an unknown pure 
Gaussian state pi„ which she wants to teleport to Bob. In a first strategy, Charlie is traced out and Alice and Bob implement 
a standard continuous variable teleportation protocol. In an alternative strategy (dashed detector), Bob is helped by Charlie 
who measures his mode and classically communicates the result n to Bob, who uses also this information for his local operation. 
Here we consider, for Charlie's measurement, first a local dichotomic measurement and then a local Gaussian measurement. 



III. DICHOTOMIC MEASUREMENT 

We first consider the case of a dichotomic measurement with measurement operators Eq, E\ = (I — Eq) 1 / 2 where 
Eq is an arbitrary Gaussian state with displacement do and CM Vo- This implies that for the outcome n = the 
conditional bipartite state p^ = P ( ^ 1 Tt c (EqpEq) is still Gaussian, while for the other outcome it is not Gaussian. 
One can prove (see Appendix A) that the CM of the reduced state p^ is 

y(p) — ytr _ ( EME T EMD t \ 

V ~ V { £)M E DMD ) (4j 
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where M is the following 2x2 "measurement matrix" 

M = -n [2(det Vo + 1/4) Vo + 4(det Vfc)C1 ft T (5) 
5 

with 

.9 = 4 det V det C + 2(det V + l/4)Tv{V VLCn T ) + (det V + 1/4) 2 > (6) 

and 

If we use Eq. Q in the right hand side of Eq. © instead of Eq. © , we obtain 

r (0) = ptr _ S T MS (g) 

where 

•£ = E T R - D T . (9) 

Now, the conditional teleportation fidelity is given by p(°) = (detr^))- 1 / 2 and satisfies the following 

Proposition 1 The conditional fidelities PW corresponding to the outcomes n = 0, 1 and the assisted fidelity F 
satisfy the inequality 

pW < p <F tr < P (0) (10) 

Proof. The proof of F tr < Fw is based on Eq. JHJ. Matrix M is real, symmetric and strictly positive and, since 
£ is real, the matrix S T MS is real, symmetric and positive. Likewise T tr and are real, symmetric and strictly 
positive, and from linear algebra it follows that det(r(°>) = det(T* r - E T MS) < det(r tr ). The proof of FW < F tr 
is a consequence of the previous result. The characteristic functions <I> tr , < &^ - ) and <I>( 1 * ) of states p tr , p' ' and p^ 
are related by = Pf 1 [$* r — P <I>( )] which also shows that is not Gaussian. On the other hand, the fidelity 
P« can be expressed in the form 7] FW = tt" 1 / d 2 X |$ in (A)| 2 [$W(A*, A)]* exp(-AJ( 1 )* + A*<5^) where $ in is 
the characteristic function of the input state and 5^ is an additional shift optimizing P^ 1 '. Now, it is easy to prove 
that, for every S (1 \ one has pW < P{ l [F tr - P P (0) ] < P^ l [F tr - P F tr ] = F tr . Finally, the effective fidelity 
is given by P = P P' ' + PiP^ and, since P'°) > fW, one has P > P' 1 ^. On the other hand, from inequality 
P (1) < Px l [F tr - P P (0) ] one can derive P < P F^ + [F tr - P P (0) ] = F tr . ■ 

Eq. Hll)[) shows that the dichotomic measurement leads, on average, to an assisted fidelity P which does not 
outperform F tr ', proving that the present dichotomic scheme does not seem to bring advantages in a real teleportation 
network. However the situation is very interesting from the point of view of the conditional teleportation fidelities 
p("). In fact Eq. (|10|l shows that teleportation fidelity always increases if Charlie performs a measurement and the 
corresponding conditional state is still Gaussian (p(°) > F tr ), while it always decreases with respect to the trace 
case for the outcome corresponding to the non-Gaussian conditional state (F^ < F tr ), even if in this latter case the 
conditional bipartite state is more pure than that without measurement |l5| . This result suggests which is the right 
kind of measurement to be considered at Charlie's site (local Gaussian measurement) and it will be the starting point 
of the next section IV. 

Moreover, the dichotomic scheme can be used in a probabilistic way, i.e., selecting only the Gaussian outcome. 
In this case Bob asks Alice to perform the Bell measurement and the classical communication only if Charlie's 
measurement has given the Gaussian outcome. In such a case the assisted fidelity F is just the conditional one F^°\ 
but the protocol has a success probability equal to Pq. 

We have also explicitly verified (see section V.B) that the Gaussian outcome n = can give F^ > 1/2 for the 
teleportation of coherent states when p tr is not entangled (and therefore F tr < 1/2). In other words, Charlie can 
conditionally generate remote bipartite entanglement between Alice and Bob if the dichotomic measurement selects 
the Gaussian outcome. All these considerations make clear why it is profitable to optimize the "Gaussian" conditional 
fidelity p(°) upon the measurement parameters and exactly such optimization work concerns the remainder of this 
section. 

Thus we restrict to the Gaussian outcome (n = 0), and look for the optimal Gaussian state Eq (i.e. the optimal 
CM Vo) which maximizes the fidelity F^°K As a first result we can prove the following 
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Proposition 2 For every Gaussian state Eq, there exists a pure Gaussian state E . p such that > F^ . 

Proof. For every Gaussian state Eq, there exists a Gaussian unitary transformation U such that Eq = Up(tit)U' 
where p{tit) is a thermal state with nx > mean number of photons 16] . Thus, we can rewrite the reduced state 
= p- 1 Tr c (E pEl) = P _1 Tr c [p(n T )L>V^p("'T)] so that the fidelity P (0) achieved from the tripartite Gaussian 
state p and the measurement operator Eq is the same which is achieved from the tripartite Gaussian state W pU 
and the measurement operator p(nx) i.e. F^ = F(p,Eo) = F[W pU , p(nT)]- On the other hand, denoting with 
M(n^) the measurement matrix corresponding to a thermal state p{nr), it is easy to prove that M(0) — M(ut) > 
Vn T > 0. From this relation and Eq. ©, one has that F[U^pU, p(n T )} < F[U^pU, p{0)], but F^ = F[U^pU, p(n T )} < 
F[U^pU,p{0)} = F[p, Up(0)&] with Up{Q)U^ = £ , P pure Gaussian state. ■ 

According to this latter result, the optimal Gaussian measurement operator Eq is actually a projection onto a pure 
Gaussian state and therefore has to be searched within the set of squeezed states |a, e) = D(a)S(e) |0). Here D(a) is 
the displacement operator with a complex amplitude, while S(e) is the squeezing operator with e = r exp(2itp) and 
r, ip are the squeezing factor and phase respectively [17]]. Since the CM of the input state, T^„, is given, Charlie has 
to optimize the protocol only with respect to the CM Vq of the squeezed state \a,e), which is given by 

, , 1 / ^shrV + C^cos 2 ^ (£- 1' 1 ) cos tp sin tp \ , 
Vo\£,<P)- 2 ^ (^_ e -i)cos^sin^ £cos 2 tp + sin 2 tp J {LL) 

where £ = exp(2r), and therefore the optimization has to be done with respect to £ and tp. Finding a global maximum 
point (£, ip) is difficult in general and therefore we split the problem in two steps: we first maximize the fidelity 
p(o,p) = F(^,tp) with respect to £ for an arbitrary but fixed ip, and then we maximize the result F(£(tp), tp) with 
respect to ip. Mathematically speaking the function F(l;,ip) is bounded and continuous in the domain ]0, +co[x [0, 7r] 
and we implicitly have to consider its continuous extension in [0, +oo] x [0, 7r] in order to surely have the existence of 
global extremal points. The first step of maximization is solved by the following 

Proposition 3 For every squeezing phase tp e [0,7r], Charlie can select a squeezing factor £(tp) such that F {^(ip) , p) > 
-F(£, tp) V£ (phase-dependent global maximum point). The point £,(tp) can be derived analytically from the CMs V and 
Vi n , according to the following four-step procedure: 

1. Construct the matrices T tr of Eq. E of Eq. and U = S^r t 'TJ T E T . 

2. Define a 2-D vector 



dctC+ 1/4 
(det S) 2 - Tr{QCn T U) 



(12) 



and a 2-D phase- dependent vector 



CM) 

where tp 1 " = (smtp, cos tp). 
3. Define the scalar product 

■y{tp) = u-k{tp) (14) 

and the third component of the vector product 

cu(p) = ^[k(p)xk(p-n/2)} z . (15) 

4- Denote with p{tp) the p-dependent logic proposition 7(95) < A j{ip — tt/2) < 0. 
Then: 



p(tp) = 1 <=> £{(p) = * — = f_ {tp) (16) 

j{tp - tt/2) 

p{tp) = <=> l{tp) = V £» = +oc. (17) 
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See Appendix B for the proof. 

The second step concerns the maximization over the squeezing phase ip. Using vectors u T = (u x ,u y ) in (|12f> and 
k(ip) T = (k x (ip),k y ((p)) in (jT3|l . the fidelity can be written as 



i , r tr ~% + (g/^Mg " 7[/2) + (TVg)^) 1 " 1 ' 

«* + (e/2)fc,(^ - tt/2) + (£-i/2)k y (v) . 



(18) 



We then consider the piecewise continuous function of ip, £ = defined according to (|16|l . I|17|) and the corresponding 
phase-dependent teleportation fidelity F(<p) = F(£((p),(p) which is continuous on [0, tt]. From Eq. (| 1 If) one has 

V (t, <p) = Vo(^ ^ + tt) = VoiC 1 ,* + t/2) (19) 

and therefore F(0, y?) = F(+oo, ip + tt/2) and .F(0, <p + tt/2) = i^+oo, ip). This implies that finding the maximum 
point (p of F(ip) is equivalent to find the maximum point of the piecewise continuous function 

* m -\F(0,ip) = [detr tr -k x (ip)/k y (cp)]- 1 / 2 if p(p)=0 . 1 Uj 

Now, one has three cases: i) <p is a stationary point of F(0, ip); ii) <p is a stationary point of -F (£_(</?), (/j); iii) </? is one 
of the border points dividing the intervals where p(<p) = 1 from those where p(ip) = 0. We report a simple analytical 
expression of the final global maximum point only in the first case, while in the other two cases the expressions are 
extremely involved. In case i), defining the 2x2 matrix r = UflCfl T , the stationary points ip± of F (0, cp) are given 
by the relation 



O T 12 ~ T 21 ± T ll ~ T 22 )y/(T 11 ~ T 22 2 + 4Tl 2 T 2 l 

cos2<p± = ■ . (21) 

{Til - T22Y + ( T 12 + 721 ) 

In many cases of practical interest (for instance when coherent states or ip = squeezed states are teleported through 
a CM V with diagonal blocks, as for example in 0, 0, H ) > the above procedure allows to find the maximum point 
te(ip),(p) and the corresponding optimal conditional fidelity Fmlx — Fte(ip),<p) quite quickly. In some easy cases 
when matrices U and C are proportional to the identity, we see from l|13|) that the above optimization becomes 
(/^-independent and therefore the maximum point is given by £ = 1 if 7 < or by £ = if 7 > 0. In the first case 
the optimal Gaussian state is a coherent state, i.e. E^ vt — \a) (a\ (with a arbitrary), while in the second case it is 
an infinitely squeezed state, i.e. E^ pt = \X(ip)) (X(ip)\ where X(tp) = 2~ 1 / 2 (ce~ lip + $ e ttp ) (phase tp and eigenvalue 
X(tp) are arbitrary). 



IV. LOCAL GAUSSIAN MEASUREMENT 



The above optimization results (propositions 2 and 3) refer to the conditional scheme where Charlie performs a 
dichotomic measurement and the Gaussian outcome n = is selected, and they can be directly used for the conditional 
generation of entanglement in that configuration. However, these results can be extended to a different scheme where 
all the measurement outcomes are Gaussian and therefore all the conditional fidelities can potentially outperform 
F tr according to Proposition 1. More in detail, Proposition 1 suggests to consider for Charlie a measurement which 
creates a conditional bipartite state which is Gaussian for every outcome, and we surely achieve this condition if 
we consider for Charlie a local Gaussian measurement, i. e., a local measurement {E(n)} transforming a Gaussian 
multipartite state into another Gaussian state for every measurement outcome n. Notice that here we include in the 
Gaussian states also asymptotic Gaussian states, such as the infinitely squeezed states. Examples of local Gaussian 
measurements are provided by heterodyne measurement on a single mode c, i.e., {\a)(a\/y/TT,a € C}, or they are 
obtained when the c mode is coupled to ancillary modes by a Gaussian unitary interaction and then the ancillas are 
subject to heterodyne or homodyne measurement. Consider then an assisted protocol where Charlie performs a local 
Gaussian measurement {E(n)} on his mode c and classically communicates the measurement result n to Bob (who, 
in turn, makes a drift-cancelling displacement depending upon the measurement outcomes of both Alice and Charlie). 
It is possible to prove a result analogous to Proposition 2 of the dichotomic case: 

Proposition 4 For every local Gaussian measurement {E(n)} with fidelity F, there exists a "pure" local Gaussian 
measurement {E e (a) = \a,s) (a,s\ /y/n,a € C} with a suitable e, such that its fidelity F(e) > F. 
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Proof. Suppose that Charlie performs an arbitrary local Gaussian measurement {E(n)} on his mode c, so that 
the conditional reduced state of Alice and Bob is given by p^ = P(n)~ 1 Tr c (E(n)pE(n) il ) corresponding to a fidelity 
(the effective fidelity of the protocol is the average over the results F = E n - P ( n ) i?( " ) )- Suppose now that 
Charlie performs a further local measurement on c given by a dichotomic measurement and the Gaussian outcome Eq 
has been selected. In such a case the reduced state will be oc Tr c (E E(n)pE(n^El) = Tr c (E( n y E^E(n)p) = 
Tr c {E{n)' pE[n)'^) where E(n)' = [E(n)^ E^E^)] 1 / 2 is a Gaussian state 0], and the corresponding fidelity will be 
From Eq. (HUJ) we have F^ n > < F^' Vn, implying F < E„ p (") i?(n) ' < where is the maximum 

value achieved for a particular outcome n. Now, Proposition 2 tells us that there exists a pure Gaussian state 
E p = \a, e) (a, e\ such that F^ > F^' , with a suitably chosen squeezing complex factor e, while a can be arbitrary. 
It is then evident that one has a teleportation fidelity F^ also if Charlie directly applies the measurement E e (a) = 
\a,e) (a,e\/y^TT on the tripartite state p and classically communicates the result a to Bob, so that one has F(e) — 
F<j>) > > F. ■ 

Trivially the previous proposition assures the existence of a local Gaussian measurement of the pure form {E e (a) = 
\a, e) (a, e\ /^/tt, a € C} which leads to an assisted fidelity F(e) > F tr . In fact it is sufficient to consider {E(n)} = I 
and apply the proposition. More importantly it implies that the optimal local Gaussian measurement must be searched 
for within the set of pure measurements {{E £ (a)}, e £ C}, which is equivalent to maximize with respect to the 2x2 
CM Vb(e) = Vb(£, (p) of Eq. (|TT)l . Thanks to this result, the optimization procedure is exactly the one given for the 
dichotomic measurement, i.e., it is given by the maximization over the squeezing factor as in Proposition 3 and by the 
subsequent maximization over the squeezing phase. Repeating such procedure it is possible to find an optimal pair of 
parameters V?) which describes the optimal local Gaussian measurement {Eg(a)} and provides the corresponding 

optimal assisted fidelity F(£, (p) using l(T%|l . Notice that for finite squeezing (£ ^ 0, +oo) the measurement {E g (a)} can 
be realized by first applying a unitary squeezing transformation 5(e) to mode c and then making heterodyne detection. 
For infinite squeezing (£ = 0, +oo) the measurement {E & (a)} is instead equivalent to a homodyne detection, i.e., to 
\X(0)) c {X(0)\ for £ = +oo and to \X((p + tt/2)) c {X(<p + vr/2)| for £ = 0, where X(tp) = 2~ x / 2 (c e'^ + gt e itp ). As 
for the dichotomic case, we can use this optimized measurement to create conditional bipartite entanglement between 
Alice and Bob and now this can be done in a deterministic way since all the outcomes are Gaussian. Note that this is 
not in contrast with the impossibility of entanglement distillation of Gaussian states with local Gaussian operations 
and classical communications |13| because here we only have a transfer of entanglement resources from a tripartite to 
a bipartite state. In the following section we give explicit examples of application of our optimization procedure. 



V. EXAMPLES 



A. Optimization of fidelity 

As an example of application of our theoretical results, we consider a three- mode Gaussian state with CM 

/ ql wR wR \ 

V = wR si tl (22) 

y wR ti si ) 

where R = diag(l,— 1), / is the 2x2 matrix identity, and the coefficients q, s,t and w are real numbers. The CM 
(12211 represents a genuine CM (i.e. it corresponds to a physical state) if and only if [3| 

V- l -J>Q (23) 

where 

j= I n \ (24) 




with f2 given in J7|). Setting s = (q + l)/2, t = q/2 and w = [yj (2q — l)(q + l)]/2 in J22J, we have a q— dependent 
CM V{q) which is genuine for every q > 1/2. This is exactly the correlation matrix considered in where a novel 
(cooperative) telecloning protocol is proved to outperform the standard (non cooperative) one for increasing values of 
q. We call q the 'noise parameter' since it determines the linear entropy of the bipartite reduced state p tr = Tr c (p) 
used by Alice and Bob in the non-assisted protocol, i.e., S(p tr ) = 1 — Tr a h(pj r ) = q/(q + 1). 
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Suppose that Alice, Bob and Charlie possess modes a, b and c respectively, and Alice wants to teleport a coherent 
state (Vi n — 1/2) to Bob with the help of Charlie. Due to the simple form of V and l^ n , it is straightforward 
to compute the optimal local Gaussian measurement which Charlie can perform in order to optimize teleportation 
fidelity. From Proposition 4 we know that it has the pure form {|a, e) (a, e\ /y^, a G C} and applying Proposition 3 
we can calculate the value e = (£, (p) corresponding to a maximum. The procedure of Proposition 3 goes as follows: 



i. r* r = hi with h 

S = (w - t)I, 
U=(w- tfhl. 



(l + q + s-2w), 



2. u 



s 2 + 1/4 
(w-t) 2 [(w-t) 2 -2sh] 

{w-tfh \ _ 
s 



= k. 



Note that, since Charlie's submatrix si and matrix U are proportional to the identity, vector k((p) becomes 
^—independent and therefore all the subsequent optimization procedure becomes ip— independent. 



3. 7 = (u; - t) 2 [s(w - t) 2 - s 2 h ■ 
w = 0. 



h/A], 



4. As already pointed out at the end of section III, we have the following simplification for Charlie's optimal 
measurement 



7 < 
7 > 



(25) 
(26) 



The first case corresponds to an heterodyne detection while the second case corresponds to an homodyne 
detection. 

We have 7 < for every q > 1/2 (see Fig. |2} and therefore heterodyne detection (£ = 1, <p arbitrary) is the optimal 
local Gaussian measurement which Charlie can perform in order to maximize teleportation fidelity of coherent states 
between Alice and Bob with this kind of shared channel. The corresponding optimal assisted fidelity F(£,(p) is 
computed from 118|) setting £ = 1 and tp arbitrary, and it is given by 



F = {h 2 — (w — t) 2 (s + l/2)- 2 [(2s + l)h -(w- t) 2 ]}- 1 ' 2 



(27) 



which is a function of the noise parameter q. Fig. [31 clearly shows the improvement provided by the optimal assisted 
fidelity F{q) with respect to the non-assisted fidelity F tr (q) = h^ 1 for every q and especially for increasing noise in 
the channel. 
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FIG. 2: </?-independent scalar product 7 versus the noise parameter q for 1/2 < q < 50. 
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FIG. 3: Optimal assisted fidelity F(q) (solid line) and non-assisted fidelity F tr (q) (dashed line) versus the noise parameter q 
for 1/2 < q < 50. 



B. Conditional generation of entanglement 

Consider now a three-mode Gaussian state p with CM 

/ al fl eR \ 

V = \ fl bl dR ] (28) 

\ eR dR cl J 

where R = diag (1, — 1), J is the 2x2 matrix identity and 

{a, 6, c, d, e, /} = {10.15, 5.52, 15.2, 8.87, 12.3, 6.96}. (29) 

Numerical values in (|29|) are chosen so that CM (|28|) represents a genuine CM. Tracing out mode c (Charlie), the 
remaining bipartite state p tr = Tr c {p) of modes a and b has CM 

V tr - ( fi Ij ) (30) 



and one can verify the separability condition 2C 



where 



V tr - l -J>0 (31) 



( " o ) <*->) 



Condition (|31|l means that the reduced state p tr shared by Alice and Bob is a separable state and therefore it cannot 
allow a quantum teleportation, i. e. it leads to a fidelity F tr < 1/2 for teleportation of coherent states. 

Here we give an explicit example where Charlie can conditionally create bipartite entanglement between Alice and 
Bob by performing an optimal Gaussian measurement at his site and then communicating the result. Consider, 
for simplicity, teleportation of coherent states, and compute the optimal local Gaussian measurement applying the 
procedure of Proposition 3: 

1. T tr = (a + b + 1)1 - 2fR, 

E = el - dR, 

(e-d) 2 (a + b+l + 2f) 

{e + d) 2 (a + b+l-2f) 



U 



, c 2 + l/4 

- " - 1 ( e 2 - d 2 ) 2 - 2cA 

A + Ai cos(2(/j) 



%>)= , 

where A = (e 2 + d 2 )(a + b + 1) - idef and Ai ee -2/(e 2 + d 2 ) + 2de(a + 6+1). 
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3. 7 (fp) = c(e 2 - d 2 ) 2 - (c 2 - 1/4) A + (c 2 + 1/4) Ai cos(2p), 
u)(tp) = cAicos(2(^). 

The next step concerns the study of the value of the proposition p(<p) : j((f) < A j((f — 7r/2) < 0. Fig.0]shows that 
7(93) < for every t/j, and therefore p(<p) = 1 for every tp, so that the phase-dependent global maximum point £(<p) is 
always given by £—{<p) as in 1|16[1 . For this reason the maximization over the squeezing phase ip is equivalent to find the 
maximum point of F(ip) = ((p), ip) (see l|20() and H18I 11. which we have plotted in Fig. EI Maximum points take the 
values (pk — kir/2, k = 0, 1, ... so that we can choose (p = 0, which gives £ = £(</?) ~ 0.087 and (p) = F((p) ~ 62%. 
In conclusion, Charlie's optimal local Gaussian measurement is equivalent to a squeezing transformation S(e) of his 
mode given by e = <^) = (0.087, 0) followed by an heterodyne detection. Such a measurement (and the subsequent 
classical communication) implies a fidelity of 62% for the teleportation of coherent states and therefore it is sufficient 
to create bipartite entanglement between Alice and Bob. 
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FIG. 4: Scalar product 7 versus the squeezing phase ip. This refers to a tripartite Gaussian state having CM 1281 with parameter 
choice iJSSJ. 
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FIG. 5: Fidelity F versus the squeezing phase (p. This refers to a tripartite Gaussian state having CM 1281 with parameter 
choice ffity. 



1. Different parameter choice 

It is instructive to study a different parameter choice leading to a more involved situation. Setting 

{a, 6, c, d, e, /} = {0.55, 0.89, 0.94, 0.74, 0.249, 0.12} (33) 

in (|28() . we have again a genuine CM and a separable reduced state p tr for Alice and Bob. However in this case the 
proposition p(tp) is true only for ip €](pi,<p2[ U ]ip3,(pn[= 1Z, where the border points ipk are given by pi — 0.339, 
p>2 = 7r/2 — 0.339, tp 3 = tt/2 + 0.339, and p4 = it — 0.339 (see Fig. EJ. Maximization over the squeezing phase is given 
by the maximization of the piecewise continuous function 

p./^_/^-(v),v) if ^e^ 

* [lf) -\F(0,p) = [detr tr -k x (p)/k y (p)r 1 / 2 if pe[0,ir[-K ^ 
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(see H2O0 and p8|l1. In Fig. we see that F(£_((p),(p) does not have maximum points inside the region 1Z, while 
F(0,ip) has two stationary points ip_ — and (p + = ir/2 which fall in [0,tv[—TZ. They are exactly the ones derived 
from lj21|) as we can easily check noting that ri2 = T21 = and m — T22 < so that cos(2ip±) = =pl. It is evident 
from Fig. [7] that y_ is a maximum point while ip + is a minimum point for -F(0, ip). In order to find the global 
maximum, we have to compare F(p>-) = .F(0, (/?_) with the right and left limits of F(<p) at the border points ipk- 
Since F(pi) = F(^) = F(£_(^ 2 ), ¥>2) = F($_(^ 3 ) s Va) = 0.514, F(0,pa) = F(0,p 8 ) = 0.446, and F(0,^_) = 0.526, 
we have that tp = (p_ = 0. In conclusion, in this different parameter choice, the optimal Gaussian measurement at 
Charlie's site is the homodyne detection \X(ir/2)) c (X(ir/2)\ which implies a fidelity of 52.6% for teleportation of 
coherent states, and therefore it is again sufficient (together with the classical communication of the result) to create 
bipartite entanglement between Alice and Bob. 
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FIG. 6: Scalar product 7(1/3) (solid line) and its displaced form 7(^ — 71-/2) (dashed line) versus the squeezing phase if>. Denoting 
with (fk (k — 1, 2, 3, 4) their intersections with the </5-axis, the proposition p(<p) : j(tp) < OA "{(tp — tt/2) < is true only inside 
the region 1Z =]tpi, y>2[U]</?3, <pi{- This refers to a tripartite Gaussian state having CM l|28|l with parameter choice 




FIG. 7: Fidelities F(£-(tp),<p) (dashed line), F(0, <p) (solid line) and F(<p) (thicker line) versus the squeezing phase ip. The 
piecewise continuous function F(ip) is discontinuous at points ip2 and tps. This refers to a tripartite Gaussian state having CM 
H28fl with parameter choice l|33fl . 
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VI. CONCLUSION 



We have studied how a two-party teleportation process (between Alice and Bob) within a three-party shared 
quantum channel can be conditioned by a local measurement and a classical communication of the third party 
(Charlie). In particular our analysis has been carried out for a shared Gaussian channel, the teleportation of pure 
Gaussian states, and two general kinds of local measurement at Charlie's site. We have first shown the case of a 
dichotomic measurement and we have proved that the non-Gaussian outcome always worsens the fidelity while the 
Gaussian outcome always improves it, even allowing the conditional generation of entanglement. Then we have shown 
how the dichotomic measurement can be designed so that the Gaussian outcome optimizes the teleportation fidelity 
and we have extended such results directly to the case of a local Gaussian measurement at Charlie's site. From the 
knowledge of the correlation matrices (the one of the shared tripartite Gaussian state and the one of the state to 
be teleported), Charlie can always determine and perform an optimal local Gaussian measurement given by a set of 
squeezed states with squeezing factor £(</?) and squeezing phase ip, maximizing the fidelity of teleportation of pure 
Gaussian states between Alice and Bob. 

It is an interesting and still open question to establish if this optimal Gaussian measurement is also the best among 
all possible measurements at Charlie's site. Proposition 1 shows that in the dichotomic case, the Gaussian outcome 
yields always a better result than the non-Gaussian one. This fact and also the fact that we are here considering 
the particular task of teleporting a one mode Gaussian state employing a tripartite Gaussian state suggest that 
this optimal Gaussian measurement can actually be the best possible measurement Charlie can do to maximize this 
specific teleportation fidelity. Notice that the recent paper of Ref. [2l| has shown that, for 1 — > 2 cloning of coherent 
states, even though the joint fidelity is maximized by a Gaussian doner, the single-copy fidelity is maximized by a 
non-Gaussian doner. However, from the point of view of teleportation, the Ref. [21] gives a support to our conjecture 
that the fidelity of teleportation, for Gaussian input and Gaussian quantum channel, is optimized by a Gaussian 
measurement. In fact, in Ref. [21], the particular optimal non-Gaussian doner which realizes F\ = 1 and F 2 — (with 
Fk the single-copy fidelity of the k th clone) , and therefore gives the optimal quantum teleportation from the coherent 
input to one clone, actually coincides with a Gaussian doner (see Fig. 1 in Ref. [21] at the extremal points (1,0) and 
(0,1)). 

We finally notice that the procedure sketched in this work can be applied to all CV teleportation networks based 
on a multipartite Gaussian state and can be used also for the conditional generation of bipartite entanglement. 
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APPENDIX A: DERIVATION OF THE CORRELATION MATRIX OF EQ. © 

We can derive the CM of Eq. using the symmetrically ordered characteristic function 17]. The measurement 
operator Eq is a state of mode c and it can be expressed as 

E = 7T" 1 J Dt( Vc )<Po(Vc)d 2 Vc (Al) 

where D c (n c ) — ex.p(rj c c^ — T]*c) is the displacement operator acting on the Hilbert space of mode c, $>o(r) c ) is the 
corresponding characteristic function, and r\ c = r/^ + irp c is a complex variable corresponding to the annihilation 
operator c 0] ■ Since Eq is a Gaussian state, we have 

$o(?7c) = exp(-7^V f} c + idfiffc) (A2) 

where Vq is the CM of the state, do € M 2 the displacement, and ff£ = (rj^, —rjj?) an M 2 vector connected to n c . In the 
same way, the total three-mode Gaussian state 



p = it 3 



c)d 2 r] a d 2 r] b d 2 ri c (A3) 



is associated to the characteristic function 

$(ff) = exp(-T] T Vfj + icFrj) (A4) 
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where V is the CM of Eq. (JTJ, d S 1R 6 is the displacement, and ff = (77^, 77^) an M 6 vector connected to (rj a , r/b, T] c ) 
as shown before. 

The conditional reduced state 



corresponds to a characteristic function <&(°>(ri a , r/b) which is (by definition) 

&°\r)a,Vb) = Tr ab [p^b a (r, a )bb(r,b) 

Putting Eqs. IjAlfl and i|A3(l into Eq. (|A5|) . and the subsequent result in Eq. (|A6|) . we obtain after some algebra 

'fin* — §* k 



$o(^)$o(K)$(?7Q,?7fc,-^ - K)exp 



d 2 M 2 K 



(A5) 



(A6) 



(A7) 



Inserting now Eqs. (|A2|) and i|A4(l into Eq. I|A7() . and adopting the ]R 4 variable ^ = (-d 1 , — l & R , k 1 , —K R ), we obtain 



$ (0) (M) = ir-zp^^irj^rjb) / exp(-<^Aff+ v T <^ 



where 



$ tr (Va,Vb) = exp 



Vb 



Vb 



•*(<? P ) 



is the characteristic function of p* 



M EE 



C+ffi C + V 



(A8) 



(A9) 



(A10) 



is a 4 x 4 matrix expressed in terms of the 2x2 sub-matrices Vo, C (Charlie's submatrix in Eq. QJ) and (defined 
in Eq. 0), 



2{E T fj a + D T rj b ) + i{d a -d c ) 
2(£ T 7f a + £ T 7f b ) + z(d -4) 



(All) 



is an R 4 vector, with E and D the off-diagonal 2x2 submatrices in Eq. and d c is the displacement of Charlie's 
reduced Gaussian state p c = Tr a b(p). Solving the integral in (|A8|) . we have 



<l> (0 \rj a ,rJb) = P Q - 1 <£ tr (rJa,rib) 



exp 



>/9 



(A12) 



where g = det M is given in Eq. ©. Inserting now Eq. (|A9|) in Eq. I|A12|I . and using <J>(°)(0, 0) = 1 (<=> Tr a6 (p(°)) 
1), we get 



<£> {0) {r} ai Tlb) = exp 



+ i 



T ( ?7a 



where 



j(0) =c jtr + ( EM(d - 4) 
V £>M(<io - d c ) 



(A13) 



(A14) 



is the displacement, with M given in Eq. JSJ), and the CM V^ ' corresponds to the expression of Eq. 
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APPENDIX B: PROOF OF PROPOSITION 3 

Using vectors u T = (u x ,u y ) in ((Hj) and k(ip) T = (k x (<p),k y (cp)) in JTBJ>, the fidelity F(£,ip) can be written as in 
(I18|l . For an arbitrary 99, we want to compute the stationary points £±(y) of F(£,tp) in the £ variable. We then 
introduce the quantity 



4g 2 , 9 (g,y) 2 d 



iV(^) = ^F(^) (Bl) 



so that they are given by iV(£, 99) = (for < £ < +oo). If the stationary points £±(<p) exist (i.e. they are finite and 
positive), they have the form 



iw - n/Z) 

with 

-(99) = c^) 2 + 7(9? - ^2)7(99) (B3) 
and 7(92) and u>((p) defined respectively in (|14fl and I] 15(1. From the product 

W*H = -fP^ (B4) 
[7(99-77/2)] 

it follows that 7(99 — 77/2)7(99) > implies the existence of only one stationary point (£+(99) and £-(99) have opposite 
sign). In particular 

7(93 - vr/2) < A 7(99) < (99) > 0, (B5) 

while 

7(99-77/2) > OA 7(99) > 0=^(99) > 0. (B6) 

Since 



7V(e±(^+£,^=7(^-V2)e 2 ±2yS(^)£ (B7) 
it is easy to prove that in (|B5(1 £—(92) is a global maximum (= £ ma x(<p)) while in (|B6|I £+(99) is a global minimum 

(= £mi n (<p)). 

Looking at (|B4ll . we then analyze the other conditions 7(99 — 77/2)7(92) < and 7(99 — 77/2)7(99) = 0. 



• If 7(99 — 77/2)7(99) < then £+(99) and £-(99) are different from zero and have the same sign. Suppose that 
£±(99) > and consider the case 7(99 — 77/2) < (^7(99) > (the proof is analogous in the other case 7(99 — 77/2) > 
OA 7(99) < 0). From 1B2|I it follows that 09(99) < 0. On the other hand, from 7(99) = u x k x (ip) + u y k y ((p) > and 
7(99—77/2) = u x k x (ip—Tr/2)+Uyky((p—ir/2) < 0, it follows that k x (ip)ky((p—Tr/2) — k y (ip)k x ((p—Tr/2) — 2u>{tp) > 
which is impossible. For this reason 7(99 — 77/2)7(99) < implies the non existence of the stationary points £±(99). 
Now, since the extended function F(£, 99) surely has a maximum and a minimum, the global extremal points 
£min(99) and £ max (<p) must be the two border points or +00. 

The last condition to be analyzed, i.e. 7(99 — 77/2)7(99) = 0, must be distinguished in more cases. 

• First suppose that 7(99 — 77/2) = 7(99) = 0. Using this condition in 118|) . it leads to F(£, 92) = 

[detr tr — k x (ip)/ky(ip)] , making the fidelity £— independent, so that we can freely choose £ m ; n (99) and£ max (99) 
at the border. 

• Next suppose that 7(99 — 77/2) /OA 7(99) = 0. In such case 7(99) = u)(ip) = — \k y {(f) /2u x Yt{ip — 77/2) 
and since k y (tp) > 0, u x > we have that u)(ip) and 7(99 — 77/2) have opposite signs. On the other hand 
7(99) = ==>■ £±(99) = [u(ip) ± |w(99)|] /-f((p — 77/2). Consider now 7(99 — 77/2) > (the proof is analogous in the 
other case 7(99 — 77/2) < 0), therefore uj(<p) < and we have £-(99) = 2uj(<p)/'y(<p — 77/2) = —k y (ip)/u x < while 
£+(99) = 0, concluding that stationary points do not exist (global extremal points at the border). 
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• Finally the case j(ip — 7r/2) = A -f(tp) 7^ can be taken back to the previous ones using periodicity ar- 
guments. Setting tp = ip — it/2 in the last logic proposition, we achieve "f(p>) = A j(tp + 7r/2) 7^ 0, from 
which we can derive both "f(ip — tt/2) — 7(1^) = and 7(1^ — 7r/2) 7^ A ^{0) = 0. In either cases we have 
(6nm(<p) = OACmax(^) = +oo)V(£ m in(<^) = +oo A ^max(^) = 0). From I© we have F(£, (p) = F{£_-\p+n/2) = 
■F(£ -1 , ip) and therefore (£ min (</?) = +°° A £ max (<p) = 0) V {^ m m(v) = A £max(^) = +00). 

In conclusion we can summarize all the cases as follows: 

7 - 7T/2) 7 (^) < <=> (6ni„(<^) = +OO A Cmax(<^) = 0) V (£ m in(^) = A £ max (<£>) = +Oo) (B8) 

>y((p - vr/2) > A 7 (^) > = Z+(v>) A = V £ ma xM = +00) (B9) 

7(<p - 7T/2) < A 7 (VJ) < <S=^ Cmax(<y5) = (^) A (Cmin(<y5) = V CminM - +Oo) (B10) 

from which we derive propositions (|16J) and Ijl7|) when we consider only the global maximum point £max(y) = 
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